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i2 . ABSTRACT 

! Diffuse gamma rays probe the highest-energy processes at the largest scales. Here 

\ we derive model-independent constraints on the hadronic contribution to the Galactic 

^ \ and extragalactic 7-ray spectra at in the energy range 50 MeV < < 10 GeV. The 

' hadronic component is dominated by emission from neutral pions, with a characteristic 

■ spectrum symmetric about m^o /2. We exploit the well-defined properties of the pion 

^T) • decay spectrum to quantify the maximum pionic fraction of the observed 7-ray intensity. 

. We find that the Galactic spectrum above 30 MeV can be at most about 50% pionic. 

, The maximum pionic contribution to the extragalactic spectrum is energy dependent; it 

fH I also depends on the redshift range over which the sources are distributed, ranging from 

C^l as low as about 20% for pions generated very recently, to as much as 90% if the pions 

O ' are generated around redshift 10. The implications of these constraints for models of 

' 7-i"ay and neutrino emission are briefly discussed. 

ci '. 

J> . Subject headings: cosmic rays - gamma rays - nuclear reactions, nucleosynthesis, abun- 

S> ! dances 

1. Introduction 

The prominence of diffuse emission in the 7-ray sky above > 50 MeV has been known since the 
earliest days of 7-ray astronomy itself (Fichtel, Kniffen, & Hartman 1973). These diffuse photons 
carry unique and direct information about some of the most energetic sites and processes in nature. 
Diffuse 7-ray observations thus provide a powerful tool both (1) to test specific models of known or 
postulated astrophysical sources, and (2) to constrain, in a model-independent way, known physical 
processes which might occur in one or more sources. We take the latter approach in this paper, 
focusing in particular on the 7-ray spectrum and the constraints it places on the contribution of 
hadronic interactions to the overall diffuse background. 
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The diffuse 7-ray sky is dominated by emission from the Galactic plane (Hunter et al. 1997), 
but the presence of emission even at the Galactic poles already suggests that an extragalactic 
component is present as well (Srcekumar et al. 1998). The spectra of these two components are 
each remarkable both for what they show and what they do not show. Namely, in neither spectrum 
is there a strong indication of hadronic interactions, which are dominated by proton collisions 
with interstellar matter, which yield 7-rays predominantly through pion production and decay: 
pp pp-K^ 77. The pionic spectrum is symmetric about a peak at m-!^/2. This feature, the 
"pion bump," is notably inconspicuous in the 7-ray data. 

As we will see in detail below, the Galactic spectrum is well-described by a simple broken power 
law, with a break at ~ 0.77 GeV. No strong pion bump is observed. Hunter et al. (1997) do note 
that there is as a ~ 2a deviation in the 60 — 70 MeV energy bin, but this region in the spectrum is 
otherwise well-fit by a smooth power law. If real, this feature is remarkably narrow. Intriguingly, 
detailed models of known Galactic processes run into difficulties explaining this spectrum (and 
its simplicity). The model of Strong, Moskalenko, & Reimer (2000) includes a sophisticated 2-D 
model of the cosmic-ray, gas, and photon fields in the Galaxy, and includes hadronic interactions, 
electron bremsstrahlung and inverse Compton scattering of starlight. However, when using only 
known cosmic ray populations and spectra, this model is unable to account for the observed 7-ray 
spectrum. The spectrum above about 1 GeV is flatter than the prediction of pionic emission, so 
other sources seem to be required as well. Proposed explanations for this "GeV excess" include 
modifications to the proton spectrum, and additional inverse Compton radiation due to an extended 
halo of cosmic ray electrons (Strong, Moskalenko, & Reimer 2000). One of the main goals of this 
paper is to quantify the portion that can be pionic. 

Information about the extragalactic component of diffuse 7-rays is more difficult to obtain, 
as one must first subtract the Galactic foreground, which is large at low-and possibly even high- 
Galactic latitudes. As we will see, the nature of the extragalactic spectrum depends on the method 
used to subtract the Galactic foreground. Different techniques have recently emerged, leading to 
different results for the shape and amplitude of the spectrum. Sreekumar et al. (1998) find a single 
power-law, while Strong, Moskalenko, &; Reimer (2003) find a smaller but "convex" spectrum. In 
either case, no pion bump is seen. 

Many astrophysical sites have been proposed to explain the extragalactic emission. These 
necessarily include "guaranteed" sources, namely, active (Stecker Sz Salamon 1996; Mukherjee Sz 
Chiang 1999) and normal (Pavlidou Sz Fields 2002) galaxies. These are the classes of objects 
which have been directly detected in nearby objects, but which would be unresolved when at large 
distances. Also, there is a growing consensus that the formation of large scale structures leads to 
shocks in the baryonic gas, and thus to particle acceleration. The resulting "cosmological cosmic 
rays" have recently become the subject of intense interest (Miniati et al. 2000; Loeb Sz Waxman 
2000; Miniati 2002; Totani & Inoue 2002; Keshet, Waxman, Loeb, Springel, & Hernquist 2003; 
Berrington & Dermer 2003; Furlanetto & Loeb 2003). These can also contribute to the diffuse 
7-background, and would have emission from both hadronic and inverse Compton processes. 
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Here we wish to find model-independent constraints on hadronic and thus pionic emission 
mechanisms. We choose to focus on this component because its detection would finally confirm ob- 
servationally the theoretical expectation that the same astrophysical acceleration processes which 
give rise to non-thermal electrons (and associated synchrotron radiation) also give rise to non- 
thermal ions. Also, we wish to exploit the well-defined nature of the pion decay spectrum which 
allows us to make a roughly model-independent comparison with observations. Finally, since the 
same hadronic processes that produce neutral pions also produced charged pions and hence neu- 
trinos, our limits will have implications for neutrino production as well. 



We will consider the Galactic and cxtragalactic emission in turn. For the Galactic spectrum, 
we adopt the EGRET data (Hunter et al. 1997) for the inner Galaxy (300° <£ < 60°, |6| < 10°). 
We find that the flux density can be well-flt by a broken power law, with index —1.52 below 0.77 
GeV, and index —2.25 above: 



This simple fit somewhat overestimates the flux in the region within about ±100 MeV of the break, 
but this region will not strongly affect our results. 

Although diffuse emission from the Galactic plane dominates the 7-ray sky, the emission is 
nonzero even at the Galactic poles, which suggests that there is an cxtragalactic component. How- 
ever, it is already clear that careful subtraction will be crucial in obtaining the cxtragalactic gamma- 
ray spectrum. Several schemes have been proposed for subtraction of the Galactic foreground. The 
basic approach of the EGRET team (Sreekumar et al. 1998) is to correlate the 7-ray sky with 
tracers of Galactic 7-ray sources. The dominant source is the hydrogen column, itself derived from 
observations of neutral H at 21 cm, H2 as traced by CO, and H II as probed by pulsar dispersion 
studies. The interstellar photon field, which is up-scattered by inverse Compton processes, is also 
estimated. Sreekumar et al. (1998) find evidence for a statistically significant isotropic component, 
with flux /(> 100 MeV) = (1.45 ± 0.05) x 10~^ photons cm~^ s~^ sr~^ and a spectrum consistent 
with a single power law of index 2.1 it 0.03: 



2. 



Data 



4bs(e) = 




eoev < 0.77 
ecev > 0.77. 



(1) 




-2.1±0.03 



(2) 



where Eq = 0.451GeV and Iq = 7.32 x IQ-^cm-^sr-^s-^GeV-^ 



Recently, Strong, Moskalenko, Sz Reimer (2003) have taken a different approach in subtracting 
the Galactic foreground, based on their sophisticated and detailed model of the spatial and energetic 
content of the Galaxy. They used the GALPROP model for cosmic ray propagation to predict 
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the Galactic component and give the new estimate of the extragalactic gamma-ray background 
(hereafter EGRB) from EGRET data. Strong, Moskalcnko, & Reimer (2003) also find evidence 
for an EGRB, but with a different spectral shape, and in general a lower amplitude than that of 
Sreekumar et al. (1998). The Strong, Moskalenko, & Reimer (2003) Galactic foreground estimates 
also includes the Strong, Moskalenko, & Reimer (2000) estimate of the pionic contribution. This 
model-based constraint will serve as an important consistency check of our model-independent 
results. We used the least square method to fit their data with a cubic logarithmic function for the 
energy range 0.05-10 GeV: 

ln(4bs-E^^) = -13.9357 - 0.03271n E + 0.1091(ln Ef + 0.0101(ln Ef (3) 

In this fit energy E is understood to be in the units of GcV and I in units of photons cm^'^ s^^ sr^^ GeV 

Indeed, the latest analysis of EGRET data done by Kcshet, Waxman, & Locb (2003) also 
implies that Galactic foreground was overestimated in previous work. They find that Galactic 
foreground in fact dominates the sky and that only an upper limit on the EGRB can be placed. 
However, Keshet, Waxman, &; Loeb (2003) analysis did not contain spectral information which is 
why it is not further investigated in this paper. The data used in this paper along with the fits are 
shown in Fig.l (Galactic component) and Fig.2 (EGRB). 



3. A Simple Model for Pionic Gamma-Rays 

The general expression for the 7-ray intensity spectrum at energy e in a particular direction is 
given by the line-of-sight integral 

m = ^ / g(e,r-) ds = ^ f r(e)n^^(r) ds (4) 

where we have ignored absorption and scattering processes which are negligible for e < 20 GcV 
(e.g., Madau & Phinncy 1996; Salamon & Stecker 1998). In eq. (4) we write the 7-ray emissivity 
(production rate per unit volume) in terms of the local hydrogen density nn and the production 
rate per H-atom (Stecker 1970; Dermer 1986, e.g.,) 

= / 2 / dEp<t>{Ep) (5) 

Note that if the shape of the cosmic ray spectrum (f){E) is the same everywhere along the line of 
sight, then /(e) = r(e)A^//, where is the hydrogen column density, and thus the shape of the 
observed 7-ray spectrum 7(e) is the same as that of the source r(e). This is the case of interest to 
us. 

The production rate F reflects the production and decay of neutral pions (with cross section 
a) due to a cosmic ray flux spectrum The shape of F(e) has well-known properties that reflect 
the symmetry of the decay photons in the pion rest frame. As described in detail by Stecker 
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(1970, 1971), the underlying isotropic nature of the rest-frame emission and the cosmic-ray beam 
is encoded in the emissivity spectrum, whose only photon energy dependence is through the lower 
limit in eq. (5). This can be written as eo(e/eo + ^o/^) which clearly has a minimum at eo = 'n-7r/2, 
and is invariant under e/eo — > eo/e; these properties guarantee that the spectrum is peaked at 
Co = 69 MeV (the pion bump) and falls off symmetrically on a logF — loge plot. 

The other key property of emissivity is found in the isobar+fireball model, which provides a 
good fit to accelerator data (Dermer 1986). Namely, at high energies e S> m-,r/2, the emissivity 
goes to the power law r(e) ~ e""*" (and thus by symmetry it goes at low energies as e+°''). This 
simple asymptotic power-law dependence is what allows us to constrain the pionic contribution of 
7-ray spectra. 

Note that the region of the spectrum immediately around the pion bump depends most sen- 
sitively on the details of the pion production cross section da{Ep, E^^) jdET^ and thus on the shape 
of the proton spectrum (j)p{E) with which it is convolved. Consequently, a detection of the pion 
bump, and its width, would not only unambiguously identify a hadronic source, but would also 
constrain the spectrum of source particles. In this case, our constraints, which are based on the 
absence of a pion bump and the asymptotic behavior of the pion spectrum, become superfluous. 
We look forward to this obsolescence, due to the eventual detection of the pion bump by GLAST 
or its successors. But until then our results remain relevant. 

A convenient semi-analytic fit to the pionic 7-ray source-function was recently presented by 
Pfrommer &: Enfilin (2003). Using Dermer 's model (Dermer 1986) for the production cross section, 
they arrive at the form: 



m = e 



n[r 



p,CR 



GeV 



Sa^ \Gev) 



2e 



+ 



2e 



(6) 



+ 0.44. The effective cross 



y - 0.14a-i-6 
g4.4-2.4aT, j mbarn. Following Dermer (1986) 



The spectral index determines the shape parameter 5. 
section app they modeled to the form app = 32 x (0.96 
we take the pion multiplicity to be ^ = 2. The cosmic ray projectile number density is np{r). This 
source function peaks at half the pion rest energy. In Dermer's model the 7-ray spectral index 
a-y is equivalent to the cosmic ray spectral index i.e. o:^. = ap (Dermer 1986). Note that in our 
limits on the dimensionless fraction of observed emission that is due to pion decay, only the energy 
dependence (i.e., the shape) of the emissivity in eq. (6) is important. 

For the case of extragalactic emission, these pionic 7-rays can come from different redshifts. 
Thus, for extragalactic origin eq. (4) becomes 



T(,) - J_ / . »H,com(^)r[(l 



Z)€,Z] 



)n{z) 



(7) 



where the dimensionless expansion rate Ti-iz) = H{z)/Hq takes the form Hi^z) = y^J7in(l + zY + 
in a flat universe. The redshift dependence of the source-function V depends on the nature of the 
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emission site (galaxies, cosmological shocks, etc.). For purposes of illustration, we will use a single- 
redshift approx n{z) = nQ5{z — z*). In this approximation different z^, amount to the shift of the 
pionic 7-ray flux in log- log plot. Thus in this simplistic view the form of the source- function would 
stay the same as in equation (6) but would be substituted with E^{1 + z) where is now 
the observed gamma-ray energy. Note that in this case, any pion bump would be redshifted, and 
thus would appear at energies < m^o/2. Thus it is clear that this feature is not apparent in the 
extragalactic spectrum, which is fiat or even convex at these energies. 

Of course, any realistic case will include contributions from a range of redshifts. However, one 
can view this distribution as an ensemble of delta functions, which will be an averaging over our 
simple cases, with a redshift-dependent weighting which scales as (1 -|- z)~^nn,com{z)'H{z)~^ (c.f. 
eq. 7). 



4. Procedure 

The main goal of this paper is to place a constraint to the maximal pionic contribution to 
diffuse gamma-ray flux based on the shape of the pionic spectrum, and fact that the pion bump is 
not observed. The way to obtain this upper limit is to see how much can we increase the pionic 
contribution by changing the parameters that it depends on so that it always stays at or below the 
observed values at all energies. The parameters that we change are the projectile and target number 
densities that enter in cosmic ray production of pions and the redshift from where we assume all 
pionic gamma rays originate. The condition of matching logarithmic slopes 

cj log Jobs (^) ^ dlogI^o{E) 

dE dE ^ ' 

of theoretical pionic gamma-ray flux and the fit to the observed gamma-ray flux guarantees that 
the ratio /^o//obs is extremized (and in fact maximized for the spectra we consider). Here I^o{E) = 
nH{r)^{E) and is given in units of GeV~^ s~^ cm~^. The energy which satisfies eq. (8) thus sets 
the values of our parameters that maximize pionic flux. 

Since the energy of pionic gamma-rays depends on the redshift as stated in the previous section, 
the slope of this theoretical flux will be the following function of observed energy E and the redshift 

^log^o + ^)/m^^>f' - {2E{1 + z)/m^o)-'- 

= — a^, 7 r- (9) 

dlogE \2Eil + z)/m^o)^'' + {2E{l + z)/m^o)-^'' 

Of course, for the Galactic spectrum we take z = 0. 

The choice of depends on the origin of cosmic rays. In the case of Galactic cosmic rays 
we will be using the classic observed-i.e., propagated-yalne = 2.75 (confirmed recently by, e.g., 
Boezio et al. 2003; Alcaraz et al. 2000; Sanuki et al. 2000). For extragalactic 7-rays, the sources 
are not known, but both blazars and shocks in cosmological structure formation have received 
considerable attention. For the case of blazars, it is not clear whether the emission is due to 
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hadronic or leptonic processes. Blazar 7-ray spectral indices have a distribution which averages 
give to a diffuse flux with index ~ 2.2 Stecker & Salamon (1996); if the emission is pionic this 
would be the proton index as well. Also, it is well known that the spectral index of cosmic rays 
accelerated in fairly strong shocks is a ~ —2 (Blandford & Eichler 1987; Jones &: Ellison 1991) 
which is expected to be the case with the cosmic rays from structure formation. Although the 
spectrum of structure-formation cosmic rays is not very well known for this purpose we will adopt 
the value a-y = 2.2, which is near the strong-shock limiting value of 2, and consistent with the 
Galactic source value (see discussion in, e.g., Fields, Olive, Casse, k, Vangioni-Flam 2001), as well 
as that of blazars. 

Now we have to match the slopes of the observed gamma-ray spectra to the slope of the 
theoretical pionic flux that was given in equation (9). This amounts to equating (9) with to the 
appropriate expressions for the spectra: eqs. (12) or (14) for extragalactic, and eq. (10) for Galactic. 
We then solve for E^{z), where we put 2; = for the Galactic case, and z = for the extragalactic 
case. 



5. Results 

5.1. Galactic Spectrum 

As described in §2, we fit the EGRET data for the Galactic spectrum with a broken power-law 
(eq. 1), and we use the emissivity for a proton spectrum Up = = 2.75. In order to set up an 
upper limit to the pionic contribution we match the low-energy index —1.52 to the slope of pionic 
7-rays; fitting to the higher energy portion of the spectrum would lead to an unobserved excess in 
the low-energy portion. The logarithmic slope of galactic spectrum is then just 

= -1.52 (10) 
dlogE ^ ^ 

We now equate this with pionic slope qiven in eq. (9) and solve for E^{z = 0). This sets 
up the maximal normalization to the pionic spectrum which is plotted in the Fig.l along with the 
observed Galactic spectrum. Also plotted is the logarithmic residual function. After integration 
over energies up to 10 GeV we can finally obtain the maximal pionic fraction of the Galactic 7-ray 
flux based on the shape of the pion decay spectrum as well as the lack of as strong detection of the 
pion bump: 

/.o,Mw(> e) = -J-^^^ (11) 

where I(> e) = I{E)dE. We find pionic fraction to be /7rO,Mw(> 30MeV) = 53% and /7rO,Mw(> 
200MeV) = 81%. While this integral constraint provides a diagnostic of the hadronic "photon 
budget," we stress that the lesson of the residual plot in Fig. 1 is that the deficit is not at all 
uniform across energies, but is very large at both high and low energies. 
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5.2. Extragalactic Spectrum 

By going through the slope-matching procedure described in the previous section we can fix 
the parameters that maximize the pionic contribution to the different extragalactic 7-ray spectra 
we consider. For the Sreekumar et al. (1998) spectrum (eq. 2), the logarithmic slope is just a 
constant 

dlogE ^ ' 

On the other hand, for (Strong, Moskalenko, & Reimer 2003, eq. 3), we have 

In Jobs ^ dln(/obs^^) .,„s 

dhiE dluE ^ ' 

= -2.0327 + 0.2182 In E + 0.0305(lnE)^ (14) 

In our simplistic picture we assume that all of the pionic 7-rays originated at a single redshift. 
Thus we go through this procedure for a set of redshifts ranging from z = up to z = 10. Figure 2 
shows our maximized pionic contribution for the two extreme redshifts, along with the fits to the 
observed 7-ray spectrum and the actual EGRET data points (Strong, Moskalenko, & Reimer 2003). 
We also present the residual, which is what it is left after pionic flux contribution is subtracted from 
the observed 7-ray spectrum. Here we see that for both EGRB spectra, the residual is large at low 
energies. However, the different shapes of the two EGRB candidate spectra lead to qualitatively 
different behavior at high energies (> 1 GeV): the residual remains substantial (> a factor of 2) 
for the Strong, Moskalenko, &; Reimer (2003) fit, suggesting the need for other component(s) to 
dominate both high and low energies. But for the Sreekumar et al. (1998) fit, the residual is 
small, and thus the pionic contribution can be dominant above 1 GeV. This difference highlights 
the current uncertainty of our knowledge of the EGRB spectrum (and even its existence, Keshet, 
Waxman, & Loeb 2003). Our analysis thus underscores the need for a secure determination of the 
Galactic foreground and the extragalactic background. 

To finally obtain the upper limit for the 7-rays that originated from pion decay, we integrate 
pionic and the observed (for both fits) fiux. Then the ratio of these energy-integrated fiuxes is the 
maximal fraction of pionic 7-rays for a given redshift. 

rlOGeV , 7- / \ 

Jeo c?e/obs(e) 

In Fig. 3 we plot this ratio as a function of redshift for three different integration ranges and for 
both Strong, Moskalenko, &; Reimer (2003) and Sreekumar et al. (1998) fits to EGRET data. Note 
that the results asymptotically approach unity. A glance at Figure 2 suggests the reason for this: 
the effect of increasing the emission redshift to "slide" the pionic spectrum leftward, toward 
lower energies. As a result, the peak and low-energy falloff are redshifted out of the fit regime, 
and the remaining high-energy power-law tail of the pionic emission then provides a good fit to the 
observations. 
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6. Discussion 

We have presented model-independent upper limits on hadronic 7-ray emission based on the 
shape of observed spectra and their lack of a pion bump. Above 100 MeV, one might expect that 
gamma-rays from 7r° decay should dominate the Galactic spectrum. However, we find that they can 
make only about 50% of total Galactic gamma-ray flux. From the shape of the residual function that 
is plotted in Fig.l we can see that due to the break at 0.77 GcV, additional gamma-ray components 
are needed both above and below the ~ 250 MeV regime at which the pionic contribution can be 
maximized. The residual gamma-rays below this scale can possibly be accounted by bremsstrahlung 
and inverse Compton scattering (Strong, Moskalenko, & Reimer 2000). However, if the pionic 
component is near its maximum, there still is the need for at one more component of gamma-rays 
above ~ 300 MeV to account for the EGRET data; alternatively, a single process may be at work, 
but with pionic emission being sub-dominant at all energies. This model-independent result is 
in good agreement with the sophisticated analysis of different models by Strong, Moskalenko, & 
Reimer (2000). We also stress again that this analysis was based on the assumption that the pion 
bump is not observed. If a future 7-ray mission such as GLAST were to identify this feature that 
would allow us to set a more definite and stronger limit on the pionic fraction of diffuse gamma-rays. 

The maximum pionic fraction of extragalactic 7-rays can be seen in Fig. 3 for different methods 
of foreground subtraction. It can go from 20 — 90%, depending on the assumed redshift of cosmic 
ray origin, but also on method used in subtracting the Galactic foreground to obtain the EGRB. 
Namely, the Strong, Moskalenko, & Reimer (2003) fit gives the fraction of 20% for cosmic rays 
that originated at the present, up to about 70% for z = 10, and in both cases there is still a factor 
of > 2 deficit of high-energy photons. On the other hand, the Sreekumar et al. (1998) fit implies 
a 20% pionic fraction for recent cosmic rays and about 90% for redshift 10 cosmic rays, with no 
significant deficit at high energies.. This large variation underscores the need for a robust procedure 
for determining the EGRB, and also emphasizes the power of a firm EGRB spectrum to constrain 
emission processes. 

Our limits can be compared to the results of specific models. Miniati (2002) finds that the 
pionic component contributes about 30% of the total emission from structure-formation cosmic 
rays, with the balance arising from electron synchrotron emission. He in turn finds that the entire 
cosmic-ray component itself can be 20 — 30% of the total (Sreekumar et al. 1998) observed 
background. Thus, the Miniati (2002) model finds that the pionic contribution is ~ 6 — 10% of 
the observed level, a fraction which would be larger if the smaller Strong, Moskalenko, k. Reimer 
(2003) background were used. The work of Keshet, Waxman, Loeb, Springel, & Hcrnquist (2003) 
neglects pionic emission entirely, arguing that the synchrotron emission should dominate. 

Of course, our extragalactic constraints reflect our simple single-redshift approximation for the 

origin of cosmic rays. More precise analysis should include some kind of averaging over the redshifts. 
Further progress first requires a detailed knowledge of the redshift evolution of sources and targets. 
On the other hand, if the "pionic bump" was observed in the EGRB spectrum by GLAST or other 



-10- 



future experiments then the position of its peak could immediately tell us something about the 
mean cosmic-ray flux. This information would then give us a better handle on the star-formation 
rate as a function of redshift. 

Finally, we note that our constraints on photons of hadronic origin also have implications for 
neutrinos produced in the same processes, a connection which has been emphasized by Waxman 
&; Bahcall (1999). This is of particular interest in the case of extragalactic neutrino emission, 
which may lead to high-energy (Ei, > 1 TeV) events observable by ICECUBE (Ahrens et al. 2004). 
Waxman & Bahcall (1999) use the energy density of ultra-high-energy (> 10-*^^ eV) cosmic rays 
to derive limits on the high-energy cosmic neutrino flux. Furthermore, they find that the EGRB 
(in the power-law form of Sreekumar et al. (1998)) implies a neutrino flux above TeV which 
violates this limit by a factor of up to ~ 100, if the EGRB is of hadronic origin. Our limits on 
the pionic fraction of the EGRB are complementary to the Waxman &; Bahcall (1999) result. Our 
constraints on the hadronic origin of the EGRB are weaker than their ~ 1% fraction, but are 
derived independently, based on the EGRB spectrum itself. 

We thank Andy Strong for valuable discussions, and we are indebted to Stan Hunter for 
kindly providing us with the EGRET Galactic spectrum. We are grateful to Vasiliki Pavlidou 
for valuable discussions, particularly concerning neutrino production. This material is based upon 
work supported by the National Science Foundation under Grant No. AST-0092939. 
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Fig. 1. — In this figure we present the maximal pionic contribution to the Galactic 7-ray spectrum. 
EGRET data points are taken from Hunter et al. (1997). The lower panels represent the residual, 
that is, log[(Fii^^)obs/(-^-E'^)7ro] = log (lobs /-^tto)- Note that the kink at 0.77 GeV is unphysical and 
just due to the overshooting of the simple broken power-law fit. 
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Fig. 2. — The maximal pionic contribution to the extragalactic 7-ray spectrum, computed by 
assuming that pionic 7-rays originated at a single redshift, namely at = and = 10. EGRET 
data points for both fits taken from Strong, Moskalenko, &; Reimer (2003). Lower panels represent 
the residual function as in Fig. 1. 
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Fig. 3. — In this figure we see the maximal fraction of pionic energy-integrated flux. It is given as 
a function of the redshift of origin for the pionic 7-rays. Fluxes were integrated from Eq up to 10 
GeV. 



